The effects of polydispersity and metastability on crystal growth kinetics 
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We investigate the effect of metastable gas-liquid (G-L) separation on crystal growth in a system 
of either monodisperse or slightly size-polydisperse square well particles, using a simulation setup 
that allows us to focus on the growth of a single crystal. Our system parameters are such that, inside 
the metastable G-L binodal, a layer of the gas phase 'coats' the crystal as it grows, consistent with 
experiment and theoretical free energy considerations. Crucially, the effect of this metastable G-L 
separation on the crystal growth rate depends qualitatively on whether the system is polydisperse. 
We measure reduced polydispersity in the crystal relative to the fluid, and propose that the required 
fractionation is dynamically facilitated by the gas layer. Our results show that polydispersity and 
metastability, both ubiquitous in soft matter, must be considered in tandem if their dynamical 
effects are to be understood. 

PACS numbers: 82.70.-y, 64.75. Gh, 64.60.-i 



I. INTRODUCTION 

Substances in the category of 'soft condensed matter' 
(colloidal suspensions, polymers, proteins etc.), as well 
as having widespread industrial and medical importance, 
are interesting in their own right because they exhibit 
physics analogous to that of simpler molecular or atomic 
systems, but on much longer time- and length-scales [l|. 
The analogy stems from the importance of thermal fluc- 
tuations in the microscopic dynamics of the constituents, 
which means that the same thermodynamic and statisti- 
cal mechanical approaches can be applied in both cases 
Nonetheless, the study and exploitation of soft mat- 
ter presents a number of challenges, two of which are the 
focus of the present work. 

Firstly, a collection of mesoscopic particles which are 
nominally the same will almost certainly be polydisperse, 
i.e. will exhibit a distribution in properties such as size, 
charge, shape or chemical makeup. This is in contrast 
to e.g. a molecular fluid of pure water in which, in a 
strict sense, every molecule is identical. Efforts by various 
workers to describe the thermodynamics of polydisperse 
systems 0-[l3] have resulted so far in a reasonably firm 
understanding of the effects of a mild degree of polydis- 
persity on the equilibrium phase diagrams of hard spheres 
and related systems. Beyond this, polydispersity remains 
a pervasive but poorly understood complicating factor in 
soft matter physics - in particular, its effects on phase or- 
derin g ki netics as actually enacted in real systems are un- 
clear [LU - tl6l |. and quantitative theoretical work thereon 
is almost non-existent [13, EH • 

Secondly, and in common with simpler systems, soft 
matter's path towards equilibrium may be influenced by 
the presence of metastable states. These are the result 
of local minima in the free energy landscape which, al- 
though they do not fully minimise the free energy of the 
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system (and therefore do not appear in the equilibrium 
phase diagram), may be encountered as an intermediate 
stage. They may be long-lived, especially if the system 
must overcome some free energy barrier in order to reach 
its global minimum. A prime example is the metastable 
gas-liquid (G-L) separation of attractive particles, which 
may subsequently progress to equilibrium crystal-fluid 
coexistence [16j . The influence of metastability on phase 
ordering is a problem of general importance, with par- 
ticular application in e.g. colloid-polymer mixtures and 
protein crystallisation |L9l - l22l |. 

In this paper, we perform simulations which focus 
on crystal growth (as distinct from nucleation) in a 
model colloidal system exhibiting the metastable G-L 
coexistence described above. By varying the interac- 
tion strength and size polydispersity, we study the ef- 
fects of metastability and polydispersity on the crystal 
growth dynamics. For our parameters, we find that the 
metastable G-L separation results in a gaseous layer coat- 
ing the growing crystal, an observation which we relate to 
theoretical free energy curves and existing experimental 
work. When polydispersity is present, this layer can act 
to speed up crystal growth, in contrast to the monodis- 
perse case in which crystal growth is slowed by the pres- 
ence of the gas layer. Interpreted as a dynamical phe- 
nomenon, this suggests that fractionation at the crystal 
interface, relying on self diffusion, is facilitated by the gas 
layer. Our results shed light on the influence of polydis- 
persity and metastability on the crystal growth process, 
demonstrating the complex way in which these factors 
can interact. 

The structure of the paper is as follows. In Section ILT1 
we describe the simulation apparatus and the parameters 
used, with reference to the equilibrium phase diagram of 
our system. In Section lLUl we outline two pertinent theo- 
retical aspects of the work, viz. the free energy landscape 
of the system, and the diffusive growth of split interfaces 
of the kind which our simulations exhibit. In Section Hvl 
we present the simulation results and discuss the effects 
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of metastability, polydispersity, and the two combined 
on the crystal growth process. Finally, in Section [Vj 
we outline and briefly test a possible explanation for the 
complex interaction between metastability and polydis- 
persity, before concluding and motivating future work in 
Section ED 



II. SIMULATION 

A. Kinetic Monte Carlo algorithm 

The simulation is built on a Kinetic Monte Carlo 
(KMC) algorithm, in which the available trial moves 
are limited to small stochastic 'hops', representing the 
inertia-free Brownian motion of colloidal particles dis- 
persed in a solvent. Therefore, in contrast to equilibrium 
MC methods, KMC is suitable for studying the dynamics 
of phase transitions since unphysical moves such as clus- 
ter rearrangements, particle resizing etc., which are help- 
ful in speeding up equilibration, are not used. In common 
with comparable Molecular or Brownian Dynamics meth- 
ods, hydrodynamic interactions are neglected. Further 
details of our implementation are available in Ref. 23 1. 



B. Parameters, crystal template 

The system studied consists of TV = 5000 spherical 
particles with diameters d n drawn from a Bates (pseudo- 
Gaussian) distribution of polydispersity (defined as the 
ratio of the standard deviation to the mean) a = 
('monodisperse'), a = 0.03 or a — 0.06 [24|. The mean 
hard core diameter (d) sets the length unit, and the time 
unit is the mean time td taken for a free particle of diam- 
eter (d) to diffuse a distance equal to its own diameter. 
The hard particle cores are surrounded by pairwise ad- 
ditive square wells of mean range A = 1.15 and depth u. 
The interaction range between particles i and j depends 
on the individual particle diameters in a scalable fashion 
MM: 



V(r) 




if dij < r < Xdij 



(1) 



with d^ = (di+dj)/2. The parent volume fraction (here- 
after also referred to as 'concentration') is <fi p = 0.34. In 
the monodisperse limit, the values of <j) p and A are such 
that the equilibrium state is a coexistence of crystal and 
gas (25[. However, by choosing u (and therefore the effec- 
tive temperature T e g — l/u) appropriately, the system's 
state point can be made to lie cither inside or outside 
the metastable gas-liquid binodal, as shown in FIG. Q] 
(reproduced from Ref. (25[). Given that the critical effec- 
tive temperature corresponds to u = 1.72, we use a well 




FIG. 1. Phase diagrams reproduced with permission from 
Ref. [H for (a) A = 1.5, (b) A = 1.25 and (c) A = 1.15 
in the density-temperature plane. With our choice of units, 
T — T c ft = l/u and p = 6<j}/ir. As A decreases, the gas- 
liquid binodal becomes metastable with respect to crystal- 
solid separation. The relevant diagram for our system is (c). 



depth of u — 1.82 inside the G-L binodal and u = 1.54 
outside, where the energy unit is fc^T. In soft matter, 
one is often free to vary the interaction strength in this 
way, for instance by adding more or less polymer into a 
colloid-polymer mixture. 

Since our aim is to study the process of crystal growth, 
not nucleation, we introduce crystallisation artificially by 
placing a 10 x 10 template of immobilised particles, ar- 
ranged in the (100) face of the FCC lattice, at one end 
of the simulation cell. The simulation cell is long in 
the x axis (L x » 66), and relatively short in the (pe- 
riodic) y and z dimensions {L y = L z « 11), so that the 
growth of the templated crystal along the x axis can be 
studied for a long period of time. The template's lat- 
tice parameter is set to match the volume fraction of 
the equilibrium crystal 25]. Increasing the y and z di- 
mensions of the lattice had no measurable effect on the 
results presented. Crystalline particles are identified by 
constructing bond order vectors q 6 from the spherical 
harmonics Y 6m [12, [2f| : Particle i is flagged crystalline 



if ' > 8-5, where the sum is over the 

iV(,(j) neighbours within -\/l.45rg of particle i, ro being 
the average separation between particle i and its closest 
6 neighbours. 

After initialising the system as an amorphous fluid, 
the square well attraction is turned on and the template 
positioned, defining t = 0. An illustrative simulation 
snapshot shortly thereafter is shown in FIG. [2] In prin- 
ciple, the nucleation of independent crystals in the bulk 
fluid is possible, but we observe no such nucleation in the 
simulations presented here. 



III. THEORY 

A. Free energy landscape 

Although this work is concerned with the dynamics 
of phase transitions, it is possible to gain a surprising 
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FIG. 2. Illustrative snapshot of a monodisperse system with 
u = 1.54 at t ~ 100, showing the crystal templating method. 
The FCC template is positioned at x = 0, and in this case 
1 or 2 further crystalline layers have so far been deposited 
from the bulk fluid. Approximately half the length of the 
simulation cell is shown. Simulation snapshots were produced 
using OVITO [13|. 

level of insight into the expected behaviour from ther- 
modynamic considerations alone [28l-l30l[. by examining 
the free energy landscape of the system. To that end, 
we now calculate approximate fluid and crystal free en- 
ergy curves for our system, in the monodisperse case, by 
perturbation to a hard-sphere reference system. 
The free energy density / is given by: 

/ « /hs + y / ^r 2 g m (r)U(r) dr (2) 

where U (r) — —u and p is the number density. The in- 
tegral is over the square well range. For the hard-sphere 
contribution /hs, we use the Carnahan-Starling free en- 
ergy densit y 1311 for the fluid, and Hall's expression for 
the crystal [32(. For the hard-sphere radial distribution 
function ,9Hs( r )j we use the Percus-Yevick expression in 
the fluid |33| and that of Choi et. al. for the crystal [34| . 

FIG. [3] shows the resulting free energy densities plot- 
ted as functions of volume fraction for the u = 1.82 case 
(inside the G-L binodals, for our value of <j> p ). Allowed 
coexistences are given by the common tangent construc- 
tion. Although quantitative agreement with the simula- 
tion data of Ref. [25[ (in terms of the coexistence vol- 
ume fractions of the gas and liquid phases) is relatively 
poor, some important qualitative features are present. 
The common tangent linking the minimum of the crystal 
branch with the gas is lower than that linking the liquid 
and gas - the crystal-gas coexistence therefore has a lower 
overall free energy and is hence the equilibrium state, 
while the gas-liquid coexistence is metastable. Further- 
more, there is no common tangent between the crystal 
and liquid, which means that on the basis of free energy 
considerations, the crystal cannot locally coexist with the 
metastable liquid. 

The lack of a crystal-liquid common tangent implies a 
growth scenario like that proposed for the experimental 
colloid-polymer mixture observations in Refs. [30|, HH, 
in which growing crystallites are coated by a layer of 
colloidal gas because of their inability to coexist with 
the metastable liquid. This 'boiled-egg crystal' should 
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FIG. 3. Free energy density / as a function of volume fraction 
<f> for a monodisperse system with u — 1.82 and A = 1.15. The 
solid curve shows the fluid branch of the free energy while the 
circles show the crystal branch. The dotted line indicates the 
common tangent for metastable gas-liquid coexistence, and 
the dashed line shows the tangent for the equilibrium crystal- 
gas coexistence. There is no common tangent between the 
liquid and crystal phases. 

deplete the surrounding bulk fluid until the required 
(equilibrium) crystal-gas coexistence is achieved overall. 
Our simulation setup allows us to model the effects of 
a metastable gas layer on the growth of a single crystal 
in a system for which we have a quantitative calculation 
of the corresponding free energy landscape, allowing new 
insight into some of the physics described in Refs. (30l . [35j j . 



B. Split interfaces 

The expected growth scenario of the 'boiled-egg crys- 
tal' described in Section IIII Al is that of a 'split inter- 
face'. That is, the crystal should form an interface with a 
gaseous region, which in turn forms an interface with the 
metastable liquid, since the crystal and metastable liq- 
uid cannot coexist on the free energy grounds discussed 
above. Such interfaces and a method for predicting their 
growth rates are described in Refs. (3(| |37|, wherein the 
example used was of a crystal-liquid-gas interface as op- 
posed to the crystal-gas-liquid interface we expect here. 
We now make use of that theory to predict the evolu- 
tion of the split interface that should be formed when 
our system is inside the G-L binodal. 

As shown in FIG. HJ we assume the growth scenario 
to be such that the phases on each side of the interfaces 
are at their 'correct' densities according to the equilib- 
rium phase diagram (FIG.[T}, the densities pa, Pb, Pc, Pd 
corresponding respectively to: the metastable liquid, the 
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metastable gas, the stable gas, and the stable crystal. 
The assumption is therefore one of local equilibrium at 
the interfaces. 

Let x\(t) be the position of the gas-liquid interface, 
and X2 (t) that of the crystal-gas interface. The respective 
fluxes at the interfaces, j\ and can then be related to 
their speeds and the densities pf. 



jl = X\{t)(pA - Pb) 

h = x 2 (t)(p D - pc) 



(3) 
(4) 



We now make the simplifying approximation that the gas 
region supports a uniform gradient between the densities 
Pb and pc, as is the case in FIG. |H and further that the 
diffusion constant therein is equal to the ideal Stokes- 
Einstein diffusion constant Dq due to the low density of 
the gas. Therefore, 



ji ~ h ~ D o 



Pb - PC 
A 



(5) 



where A = x±(t) — X2(t). Substituting into Equa- 
tions |31 and U and separating variables, with the necessary 
constants of integration given by the initial conditions 
xi(t) = x 2 (t) = 0, we find expressions for the interface 
positions: 



x x (t) = j2D a ^-t+y/2D jt 



x 2 (t) 
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(6) 
(7) 



where we have defined j3 = (ps — Pc)/(pd — Pc) and 

7 = (PB - PC')/ {PA - PB) - (PB - PC)/{PD - PC)- 

To make contact with our simulations, we read off the 
values of pi from FIG.Q]at the appropriate effective tem- 
perature T c ff = 1/1.82 K> 0.55 for the systems inside the 
G-L binodal, and use the Stokes-Einstein diffusion con- 
stant defined by our choice of time unit, Do = 1/6. The 
resulting interface positions through time are shown in 
FIG. [U and in the next section are compared with simu- 
lation data. The characteristic t 1 ^ 2 diffusive growth is ap- 
parent, with the gas-liquid interface leading the crystal- 
gas interface as expected. 

We note at this stage that the crystal-gas-liquid split 
interface scenario is very closely analogous to the crystal- 
liquid-gas interfaces described in Refs. [3a, l37[- The only 
difference is that in the present work, the 'pivot' phase 
which coexists locally with both others (and therefore 
coats the crystal) is the gas, i.e. is the equilibrium phase, 
whereas in the aforementioned work it is the metastable, 
i.e. nonequilibrium liquid phase. 



PA 



PB 

pc 



X2 



XI growth axis 



FIG. 4. Schematic density profile of the expected crystal-gas- 
liquid split interface inside the G-L binodal, where the crystal 
cannot coexist with the metastable liquid. The phases from 
left to right are the crystal, gas, and metastable liquid. The 
gradient in the gas region is drawn uniform, as is assumed in 
our calculations. 
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FIG. 5. The predicted evolution of the crystal-gas (solid line) 
and gas-liquid (dashed line) interfaces for a split crystal-gas- 
liquid interface when our system is inside the G-L binodal 
(u = 1.82, A = 1.15). 



IV. RESULTS AND DISCUSSION 

In this section, we examine the simulation results in 
terms of their time-dependent concentration profiles and 
in terms of their crystal growth rates, elucidating and 
discussing the separate and combined effects of metasta- 
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bility and polydispersity on the growth process. Then, 
in Section [V] we outline and briefly test a possible expla- 
nation for our findings in terms of fractionation at the 
crystal interface. 

We present results for polydispersities of a = 
(monodisperse), a = 0.03 and a = 0.06, with square well 
depths of u — 1.82 (inside the G-L binodal, referred to 
as 'gas-mediated') oy u — 1.54 (outside the G-L binodal, 
referred to as 'gas-free'), giving 6 state points in total. 
At each state point, 6 independent initial configurations 
were used. 



A. Concentration profiles — monodisperse 

The clearest way to visualise crystal growth over the 
whole simulation time is by plotting the time-dependent 
concentration profile. To achieve this, we plot the 'area 
fraction' aroa of particles intersecting planes at a given 
x coordinate, corresponding to a local volume fraction in 
an infinitely narrow interval Sx. Since the simulation ge- 
ometry is such that the crystal interface advances along 
the x axis, these plots are a powerful way of observing 
the time evolution of concentration differences along the 
direction of growth while averaging over the axes perpen- 
dicular to it. 

In this section, we qualitatively compare representative 
concentration profiles in the monodisperse and a — 0.06 
cases. Except where noted, the phenomena described 
occurred similarly in all trajectories of the state point in 
question. 

Let us consider the monodisperse case first. FIGs. [6] 
and [7] show example trajectories for simulations at the 
gas-free and gas-mediated state points respectively. In 
the gas-free case, the growth scenario is, as expected, 
relatively simple. The bulk fluid remains homogeneous, 
since we are above the critical T e s for G-L separation. 
The crystal template causes the growth of a single crys- 
talline region with an average volume fraction of <j) ~ 0.6, 
which is the expected equilibrium volume fraction for the 
crystal (25|. The crystal, being at a higher volume frac- 
tion than the bulk fluid, depletes its surroundings of par- 
ticles, resulting in a concentration gradient between the 
local fluid and the bulk; it is this gradient which trans- 
ports particles toward the crystal so that it can continue 
to grow. 

Inside the G-L binodal (FIG. [7]), the scenario is quite 
different. Firstly, there is clearly some G-L separation 
taking place in the bulk. Focusing next on the crystal 
template at x = 0, we can see that a region of low den- 
sity gas forms in front of the crystal almost immediately, 
shielding it from the liquid with which it cannot locally 
coexist according to the free energy considerations in Sec- 
tion llll Al This is in contrast to FIG. [51 in which the fluid 
immediately next to the crystal template retains its den- 
sity for quite some time until the growing crystal depletes 
it of particles. Note also that whereas the depleted re- 
gion in FIG. [5] fades smoothly into the higher density 



Monodisperse gas-free 
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FIG. 6. Time-dependent concentration profile along the x 
axis for one of the monodisperse gas-free simulations. The 
greyscale indicates the local volume fraction, ranging from 
(f> < 0.06 (white) to (f> > 0.676 (black). As the simulation 
progresses, the crystal advances along the x axis, depleting 
the fluid in front of the interface. 



Monodisperse gas-mediated 




x 



FIG. 7. As FIG. [6l for a monodisperse gas-mediated simula- 
tion. Gas-liquid separation takes place in the bulk, and the 
advancing crystal is coated by a distinct gas layer which ad- 
vances ahead of it. The theoretical crystal-gas and gas-liquid 
interface curves from FIG. [5] are rotated and superimposed, 
showing approximate quantitative agreement with the data. 



bulk fluid, the gas next to the crystal in FIG. [Jj forms 
a substantially sharper interface with the fluid next to 
it, indicating a distinct phase boundary between it and 
the metastable liquid. The formation of this well-defined 
gas layer shielding the crystal is consistent with the free 
energy considerations above, and with the experimental 
observations of 30]. 

The speed of advancement of the crystal-gas and gas- 
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FIG. 8. Snapshot of the crystal interface for one particular 
monodisperse gas-mediated trajectory in which the crystal 
happened to exhibit dendritic growth around an impinging 
gas bubble. Such dendritic growth was not generally ob- 
served. 

liquid interfaces is in approximate quantitative agree- 
ment with the theoretical prediction in Section MI Bl 
(FIG. [5]), which we have rotated and superimposed on 
FIG. [7J for comparison. This agreement is particularly 
satisfying given that the input densities pi required for 
the theoretical calculation in Section IIIIBI were taken 
from the equilibrium phase diagram, so that there are no 
free fitting parameters. The accuracy of the prediction 
thus gives credence to the assumptions of local equili- 
bration of the interfaces and of essentially ideal gas-like 
Stokes-Einstcin diffusion in the gas region. 

Finally, it is interesting to note that for some trajec- 
tories, the crystal can grow incomplete layers, as shown 
in FIG. [HI This is the result of dendritic growth of the 
crystal around an impinging gas bubble. Such growth 
was not found to be generally present in the other tra- 
jectories at this state point - the crystal typically formed 
complete layers covered by a layer of gas spanning the 
whole width of the simulation cell. 



B. Concentration profiles — polydisperse 

We now discuss the polydisperse a — 0.06 case. FIG. 
shows the gas-free scenario. On the simulated timescale, 
essentially no crystal growth is seen to take place, except 
perhaps for a slight 'wetting' of the template by ordered 
particles. As in the monodisperse case, the bulk fluid 
remains homogeneous, indicating that the polydispersity 
has not altered the location of the critical T e g in such a 
way as might bring this state point inside the G-L bin- 
odal. 

Moving to the gas-mediated case inside the binodal 
fFIG. ITD|) . we can see the crystal is able to grow substan- 
tially despite the polydispersity. As in the corresponding 
monodisperse case (FIG. [7]), the bulk fluid separates into 
gas and liquid regions, and a gas layer covers the crystal 
as it grows. 

To summarise the qualitative observations in this sec- 
tion: the crystal template successfully induces a phys- 
ically realistic growth process, in which the crystal de- 
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FIG. 9. As FIG. [5J for a gas-free system of polydispersity 
a = 0.06. As in the monodisperse case, the bulk fluid remains 
homogeneous. However, the crystal template does not induce 
any growth on the timescale simulated. 



Polydispersity 0.06 gas-mediated 




FIG. 10. As FIG. [6] for a gas-mediated system of polydis- 
persity a — 0.06. The bulk fluid shows G-L separation, and 
the gas layer is present in front of the crystal. In contrast to 
FIG. [9] the crystal is able to grow, although more slowly than 
in the corresponding monodisperse case (FIG. [7}. 

pletes its surroundings and forms ordered, high density 
layers. When the system is inside the G-L binodal, a gas 
layer forms in front of the crystal, shielding the crystal 
from the metastable liquid. This is exactly the scenario 
predicted for our system in Section 1111 Al and proposed 
in Refs. [30j,|35[. The presence of polydispersity a = 0.06 
substantially slowed the crystal growth overall. In the 
gas-free case the crystal showed barely any growth on 
the simulated timescale, whereas the gas-mediated case 
showed substantial crystal growth despite the polydisper- 
sity. 
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FIG. 11. (colour online) Crystallinity through time for the 
6 state points studied, in each case averaged over 6 indepen- 
dent realisations. Filled symbols: gas-mediated; open sym- 
bols: gas-free. Triangles: monodisperse; circles: a = 0.03; 
squares: a — 0.06. The standard errors are approximately 
the size of the symbols. The filled triangles (monodisperse 
gas-mediated) are shown in grey (red online) to distinguish 
them from the filled circles (cr = 0.03 gas-mediated). 



C. Monodisperse crystal growth rate — the effect of 
metastability 

Having summarised the qualitative features of the 
monodisperse and a = 0.06 simulations, we next quanti- 
tatively examine the effect of those features on the crys- 
tal growth rate. This is done by plotting the crystallinity 
(the number of particles flagged as crystalline as a pro- 
portion of the total N) through time. FIG. [TT] shows the 
crystal growth at each of the 6 state points studied, each 
averaged over the 6 independent realisations of the state 
point. 

In the monodisperse case, the effect of metastability 
alone is apparent. The gas layer in front of the crystal 
significantly slows crystal growth compared to the gas- 
free case. The likely explanation for this is that the gas 
layer inhibits transport of particles into the crystal, as 
we now discuss. 

Since the crystal exists at a higher density than the 
rest of the system, as soon as the immediate surround- 
ings are depleted of particles, particles must be trans- 
ported toward the crystal from the bulk in order for it to 
grow. This transport takes place via collective diffusion 
down a concentration gradient between the low concen- 
tration near the crystal and the relatively higher concen- 
tration away from the crystal. The diffusion is described 
by Fick's law, 

J = -DV<f> (8) 

in which J is the concentration flux, <fi is the local concen- 
tration and D is the collective diffusion coefficient. Note 



that the collective diffusion coefficient, in contrast to the 
self diffusion coefficient, increases with tfi [38| . 

We now consider the phase diagram for our system, 
as shown in FIG. [T] (c). In the absence of the gas layer, 
and assuming that the interface is locally equilibrated, 
a long-range, relatively steep concentration gradient can 
exist between cf> w 0.05 at the interface and 4> p — 0.34 
in the bulk. However, inside the G-L binodal, when the 
crystal is covered by a distinct gas layer, the relevant 
concentration gradient is that which exists across the gas 
layer itself - the gas-liquid interface beyond truncates 
the concentration gradient. The gradient in the gas layer 
is between the equilibrium gas (</> « 0.01) and the slightly 
higher concentration metastable gas at <j> » 0.05, i.e. it 
is much weaker compared to the concentration gradient 
possible in the gas-free case. 

Therefore, the low-concentration gas layer reduces 
the efficiency of particle transport toward the crystal 
because: (a) the collective diffusion coefficient D in- 
creases with concentration; (b) the gas's low concentra- 
tion means it can only support a relatively weak con- 
centration gradient in comparison to the gas-free case, 
in which a strong, smooth and long-range concentration 
gradient can exist between the interface and the homoge- 
neous bulk. Although, on thermodynamic grounds, the 
driving force for crystallisation is higher inside the G- 
L binodal, the growth is slowed due to the gas layer's 
detrimental effect on particle transport to the interface. 

The effect of the gas layer in the monodisperse case is 
comparable to that of the liquid layer in Refs. [36l . [37| . 
which was also found to inhibit crystal growth. There, 
as here, the crystal is coated by a phase which advances 
ahead of the crystal and slows down the growth of the 
crystal. It is interesting that this comparison holds, given 
that in our case the crystal is being coated by the equilib- 
rium gas, whereas Refs. [36], [37} concern a crystal coated 
by the nonequilibrium liquid. 



D. Polydisperse crystal growth rates 

We consider first the highest polydispersity studied 
here, a = 0.06. As shown in FIG. 1111 and in agreement 
with the observations in Section IIV Al the gas- free state 
point now shows only very slight growth on the timescale 
simulated, whereas the gas-mediated state point shows 
significant crystal growth (although less than in either of 
the monodisperse systems). This is a qualitative differ- 
ence compared to the monodisperse case, in which the 
gas layer instead strongly slows the crystal growth. That 
is, the presence of polydispersity qualitatively alters the 
effect of the metastable gas layer on crystal growth. 

In the a — 0.03 case, the gas-mediated system is ini- 
tially slightly faster, before being overtaken by the gas- 
free system by around t = 1000. This is suggestive 
of some kind of crossover phenomenon between what- 
ever factors are dominant in the fully monodisperse and 
cr = 0.06 cases. 
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V. CRYSTALLISATION MECHANISM 
A. Fractionation 

The results presented so far show that metastable gas- 
liquid separation can affect the rate of crystal growth, be- 
cause it results in a gaseous layer coating the crystal as it 
grows. However, the qualitative nature of that effect de- 
pends strongly upon polydispersity. In the monodisperse 
case, the resultant gas layer impedes crystal growth by 
reducing the efficiency of particle transport to the crys- 
tal. With a = 0.03, the slowing effect of the metastable 
separation is reduced, while at a = 0.06, the metastable 
separation instead strongly enhances crystal growth. It is 
this effect of polydispersity on the crystal growth mech- 
anism that we now discuss. 

We propose a possible explanation in terms of a lo- 
cal fractionation process at the interface of the polydis- 
perse crystal. In size polydisperse systems, it is well 
known that phase separation is typically associated with 
some degree of thermodynamically-driven fractionation 
in which one phase e.g. contains on average larger or 
smaller particles than another, changes its overall poly- 
dispersity, etc. Such fractionation has been predicted 
theoretically [1J3, and observed experimentally [39| and 
in simulation [23]. Of course, this kind of fractiona- 
tion pertains to the equilibrium phase composition and 
may not be fully achieved in real systems on accessible 
timescales due to kinetic factors. This is a particularly 
important consideration for the crystal phase, where par- 
ticles are essentially stationary once incorporated and, 
thereafter, the long-range particle transport required for 
fractionation is facilitated only by the presence of defects. 
In general, phase composition is expected to relax slowly 
in comparison with overall density so that fractionation 
may lag behind phase separation somewhat [17] . leading 
in the crystal phase to the 'freezing-in' of a nonequilib- 
rium composition [l8j . (Note, however, that recent work 
by us has shown that some fractionation is pos sible in 
the very early stages of gas-liquid separation [23]].) 

Nevertheless, we would expect that fractionation at the 
interface of the growing crystal, where particles are still 
mobile and may be easily exchanged with the fluid, is 
quite feasible. Indeed, the slowing effect of polydisper- 
sity on crystallisation is generally taken as evidence that 
fractionation is involved to some extent [TH, [H| . We pro- 
pose that such interfacial fractionation takes place not 
just in nucleation but during the subsequent growth of 
the crystal, so that it is facilitated by the low-density gas 
layer that forms in front of the crystal when inside the 
G-L binodal. In contrast to the collective diffusion dis- 
cussed in Section IIV CI self diffusion is hindered by par- 
ticle density, so that fractionation would be frustrated 
outside the G-L binodal, where the fluid side of the in- 
terface is at a relatively high density compared to that 
in the gas-mediated case. 

From the outset, we stress that a fulsome explanation 
of the results in the previous section requires further work 



to elucidate and test other possible contributing factors, 
such as the influence of polydispersity on the equilibrium 
phase diagram. The proposals here are motivated by the 
dynamical nature of our simulation approach and by the 
intuition that the increased diffusivity of particles in the 
gaseous layer would help any dynamical sorting process 
at the crystal interface, a consideration which we expect 
to be a significant aspect of any full explanation of the 
results. In the following, we perform further analysis on 
the simulation data in an attempt to detect the presence 
of fractionation. 



B. Polydispersity of the crystal 

The results we now present concern the gas-mediated 
and gas-free a = 0.06 simulations. To improve sampling 
in the crystal, we allowed these to run up to t ss 9000. 
By this time, the crystal interface in the gas-mediated 
simulations reached x ~ 8, while the gas- free simulations 
had only grown to x ~ 3.5. We check for fractionation by 
measuring the mean diameter (d) cr yst and polydispersity 
fcryst in the crystal. To within statistical error, 

(^)cryst 

was equal to the parental mean, (d) . 

However, changes in cr cryst were detected. In FIG. 
the evolution of cr crys t in the gas-mediated and gas-free 
crystals is shown, as a function of time and of n crys t, 
the mean number of particles in the crystal. In the gas- 
mediated case, there is a clear reduction in the poly- 
dispersity of the crystal as compared to the parental 
a = 0.06, demonstrating that the crystal is selecting a 
narrower subset of the parent size distribution, eventually 
attaining a value of er crys t ~ 0.054 at the end of the sim- 
ulated time. This is qualitatively consistent with Fasolo 
and Sollich's equilibrium calculations on hard spheres, 
in which even a small parental polydispersity was found 
to lead to reduced polydispersity in a crystal coexisting 
with a fluid. This measurement in itself is significant: 
to our knowledge, such fractionation has not previously 
been measured in the dynamic crystal growth of poly- 
disperse spherical colloids (an analogous example for col- 
loidal platelets exists in Ref. [40]), although it is proposed 
as an explanation of slow crystallisation in such systems 

MM- 

The data for the gas-free crystal is subject to signifi- 
cantly larger error, because very little crystal has actually 
formed, but appears to be consistent with a small re- 
duction in polydispersity corresponding to the frustrated 
crystal growth. In any case, it is clear that the success- 
ful growth of a crystal involves a measurable reduction 
in polydispersity from the bulk fluid, indicating a frac- 
tionation process which would be enhanced by the gas 
layer. 
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FIG. 12. (colour online) Crystal polydispersity CT cryst as a 
function of t (upper pane) and n cr yst (lower pane) for the the 
gas-mediated (black) and gas-free (grey, blue online, larger 
error bars) crystals. The error on the n cryst axis is approxi- 
mately 30. 



C. Summary 

The results above show that the crystal phase is able 
to reduce its polydispersity relative to the fluid even as 
it is growing. We have proposed that the necessary frac- 
tionation must take place at the interface and therefore, 
relying on self diffusion near the interface, is enhanced 
by the presence of the gas layer. If true, this mecha- 
nism would explain the very slow crystallisation in the 
(7 = 0.06 gas-free case, where no such layer exists. 

The a = 0.03 simulations provide an interesting in- 
termediate case, in which the crystal growth in the gas- 
mediated case is initially faster before being overtaken by 
the gas-free case. In light of our hypothesis, this may be 
because the the gas-free growth is relatively slower until 
the interfacial fluid is sufficiently depleted to allow ade- 
quate diffusion near the interface. In the gas-mediated 
growth, a gas layer is formed immediately, so that growth 
is initially faster. However, once the gas-free system has 
sufficiently depleted the interfacial fluid, the interfacial 
diffusion is enhanced and it now grows the fastest, the 
gas-mediated system lagging behind for the reasons dis- 
cussed in Section IIV CI In contrast, for a = 0.06 the 
gas-free growth is so slow (presumably due to the greater 
fractionation required) that the interfacial fluid is not 
depleted enough for such a crossover to happen on the 
simulated timescale. 

We stress that the ideas outlined in this section con- 
stitute only one potential explanation of the results we 
observed in Section HVl Much further work will be re- 
quired to determine other factors (e.g. the influence of 
polydispersity-induced changes to the equilibrium phase 
diagram), and to more rigorously test the dynamical ex- 



planation outlined here (which we expect will remain an 
important part of any more comprehensive explanation). 
Nonetheless it is clear that, whatever the mechanisms, 
the influence of metastability on crystal growth kinetics 
is far from trivial, and that its role can be qualitatively 
switched by the presence of a very mild degree of poly- 
dispersity. 



VI. CONCLUSIONS 

We have studied the kinetics of crystal growth using 
a model system in which a single crystal grows from an 
amorphous fluid. By varying the interaction strength and 
the polydispersity, we were able to investigate the sepa- 
rate and combined effects of metastable gas-liquid sepa- 
ration and polydispersity on the crystal growth process. 
Taking advantage of the pre-determined growth direction 
of the crystal, we have used one-dimensional concentra- 
tion profiles to clearly display the phase ordering scenar- 
ios when the system is inside or outside the metastable 
G-L binodal. Inside the binodal, we observed the forma- 
tion of a gaseous layer coating the crystal, in agreement 
with approximate free energy curves and with previous 
experimental observations of three-phase ordering. The 
diffusive growth of the crystal-gas-liquid 'split interface' 
was modelled theoretically, and we found approximate 
quantitative agreement between theoretical and simula- 
tion interface growth profiles. This agreement validated 
our use of two approximations in the theory, namely that 
the interfaces are locally equilibrated to the densities 
given by the phase diagram (FIG. [TJ and that particle 
transport through the gas layer is approximately ideal 
gas- like due to the gas's low density. 

The dynamic influence of this gaseous layer, in the 
monodisperse case, was to impede growth by slowing par- 
ticle transport to the crystal interface, suggesting that, 
while metastable G-L separation can enhance crystal nu- 
cleation [19] , l4lj . it may (in the ideally monodisperse 
case at least) have the opposite effect on the subsequent 
growth of the crystal if the free energy landscape is such 
that the metastable liquid cannot coexist with the crys- 
tal, as was the case here. However, introducing a small 
amount of polydispersity, corresponding roughly to that 
found in a 'near-monodisperse' colloidal system, we found 
that the gaseous layer instead strongly enhanced growth. 

The crystals grown show an overall reduction in poly- 
dispersity versus the parent, consistent with existing 
equilibrium theory and equilibrium simulations. This 
fractionation is often proposed as an explanation for slow 
crystallisation in polydisperse systems, but to our knowl- 
edge has not previously been measured in the dynamic 
crystal growth of polydisperse spherical colloids. We have 
postulated that this fractionation process, requiring self 
diffusion, is facilitated by the gas layer, perhaps explain- 
ing why this layer enhances crystal growth in the poly- 
disperse case, rather than impeding it as in the monodis- 
perse case. 



10 



It must be stressed that our proposed explanation re- 
quires further investigation and probably does not con- 
stitute the full story. Other factors, such as the effect of 
polydispersity on the equilibrium phase diagram of the 
system, must be considered - further work will be re- 
quired to establish what effect those factors have, and to 
rigorously test the explanation we outlined here. Also, 
as with other comparable simulations, hydrodynamic in- 
teractions (HI) have been neglected; more complex and 
expensive simulations would be required to quantify pos- 
sible effects of HI on the dynamics observed here. 

We have presented a detailed dynamical study of crys- 
tal growth scenarios in the presence of metastability and 
polydispersity, and shown how our findings relate to the 
free energy landscape and equilibrium phase diagram of 



the system. Whatever emerges from future work, the re- 
sults demonstrate the importance of both metastability 
and polydispersity for soft matter phase transition kinet- 
ics and, moreover, that these two factors interact in a 
complex and previously unknown manner, the causes of 
which remain to be investigated further. 
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